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Abstract 

The main purpose of this article is to provide a guide to theorems 
on global properties of solutions to the Einstein- Vlasov system. This 
system couples Einstein's equations to a kinetic matter model. Kinetic 
theory has been an important field of research during several decades 
in which the main focus has been on nonrelativistic and special rela- 
tivistic physics, i.e. to model the dynamics of neutral gases, plasmas, 
and Newtonian self-gravitating systems. In 1990, Rendall and Rein 
initiated a mathematical study of the Einstein- Vlasov system. Since 
then many theorems on global properties of solutions to this system 
have been established. The Vlasov equation describes matter phe- 
nomenologically and it should be stressed that most of the theorems 
presented in this article are not presently known for other such matter 
models (i.e. fluid models). This paper gives introductions to kinetic 
theory in non-curved spacetimes and then the Einstein- Vlasov system 
is introduced. We believe that a good understanding of kinetic the- 
ory in non-curved spacetimes is fundamental to good comprehension 
of kinetic theory in general relativity. 
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1 Introduction to kinetic theory 



In general relativity, kinetic theory has been used relatively sparsely to model 
phenomenological matter in comparison to fluid models. From a mathemat- 
ical point of view, however, there are fundamental advantages to using a 
kinetic description. In non-curved spacetimes kinetic theory has been stud- 
ied intensively as a mathematical subject during several decades and it has 
also played an important role from an engineering point of view. In the 
first part of this introduction, we will review kinetic theory in non-curved 
spacetimes and we will consider mainly the special relativistic case, but 
mathematical results in the nonrelativistic case will also be discussed. We 
believe that a good understanding of kinetic theory in non-curved space- 
times is fundamental to good comprehension of kinetic theory in general 
relativity. Moreover, it is often the case that mathematical methods used to 
treat the Einstein- Vlasov system are carried over from methods developed 
in the special relativistic or nonrelativistic case. 

The purpose of kinetic theory is to model the time evolution of a col- 
lection of particles. The particles may be entirely different objects depend- 
ing on the physical situation. For instance, the particles are atoms and 
molecules in a neutral gas or electrons and ions in a plasma. In stellar dy- 
namics the particles are stars and in a cosmological case they are galaxies or 
even clusters of galaxies. Mathematical models of particle systems are most 
frequently described by kinetic or fluid equations. A characteristic feature 
of kinetic theory is that its models are statistical and the particle systems 
are described by distribution functions / = f(t,x,p), which represent the 
density of particles with given space-time position (t, x) G 1 x R 3 and mo- 
mentum p £ R 3 . A distribution function contains a wealth of information, 
and macroscopical quantities are easily calculated from this function. In a 
fluid model the quantities that describe the system do not depend on the 
momentum p but only on the spacetime point (t,x). A choice of model is 
usually made with regard to the physical properties of interest for the sys- 
tem or with regard to numerical considerations. It should be mentioned that 
a fluid model that is too naive may give rise to shell-crossing singularities, 
which are unphysical. In a kinetic description such phenomena are ruled 
out. 

The time evolution of the system is determined by the interactions be- 
tween the particles which depend on the physical situation. For instance, 
the driving mechanism for the time evolution of a neutral gas is the colli- 
sion between particles (the relativistic Boltzmann equation). For a plasma 
the interaction is through the electric charges (the Vlasov-Maxwell system) 
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and in the stellar and cosmological cases the interaction is gravitational (the 
Einstein- Vlasov system). Of course, combinations of interaction processes 
are also considered but in many situations one of them is strongly dominat- 
ing and the weaker processes are neglected. 



1.1 The relativistic Boltzmann equation 

Consider a collection of neutral particles in Minkowski spacetime. Let the 
signature of the metric be (—,+,+,+), let all the particles have rest mass 
m = 1, and normalize the speed of light, c, to one. The four-momentum of 
a particle is denoted by p a , a = 0,1, 2, 3. Since all particles have equal rest 
mass, the four-momentum for each particle is restricted to the mass shell, 
P a Pa = —tn 2 = — 1. Thus, by denoting the three- momentum by p G M 3 , p a 
may be written p a = (—p°,p), where \p\ is the usual Euclidean length of p 
and p° = \J\ + |p| 2 is the energy of a particle with three-momentum p. The 
relativistic velocity of a particle with momentum p is denoted by p and is 
given by 

Note that \p\ < 1 = c. The relativistic Boltzmann equation models the 
space-time behaviour of the one-particle distribution function / = f(t,x,p), 
and it has the form 

(dt + ^-V x )f = Q(f,f), (2) 
where the relativistic collision operator Q(f,g) is defined by 

Q(f,g) = / k(p,q,uj)x 
Jm 3 Js 2 

x [fip + a(p, Q, u)u)g(q - a(p, q, - f(p)g{q))dpduj. (3) 

(Note that g = f in P|))- Here dw is the element of surface area on S 2 and 
k(p, q, u>) is the scattering kernel, which depends on the scattering cross- 
section in the interaction process. See [15] for a discussion about the scat- 
tering kernel. The function a(p,q,Lo) results from the collision mechanics. 
If two particles, with momentum p and q respectively, collide elastically (no 
energy loss) with scattering angle w6§ 2 , their momenta will change, p — > p' 
and q — > q' '. The relation between p, q and p' , q' is 



where 



p =p + a(p,q,u)u, q' = q - a(p,q,uj)uj, (4) 
2( ff + q°)AV-(g-ff)) 

a( ^'" )= (pO + gO)2_( w . (p + , ) )2- (5) 
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This relation is a consequence of four-momentum conservation, 

p a + q a = p a ' + q a \ 

or equivalently 

p + q = p' + q'. (7) 

These are the conservation equations for relativistic particle dynamics. In 
the classical case these equations read 

W 2 + M a = bf + kf, (8) 

p + g = p' + g'. (9) 

The function a(p,q,u) is the distance between p and p' (g and q') and the 
analogue function in the Newtonian case has the form 

ad(P,<l,u) =u- (q-p). (10) 

By inserting a c i in place of a in (j3J) we obtain the classical Boltzmann collision 
operator (disregarding the scattering kernel, which is also different). In |19| 
classical solutions to the relativistic Boltzmann equations are studied as 
c — ► oo, and it is proved that the limit as ctoco of these solutions satisfy the 
classical Boltzmann equation. 

The main result concerning the existence of solutions to the classical 
Boltzmann equation is a theorem by DiPerna and Lions [Mj that proves 
existence, but not uniqueness, of renormalized solutions (i.e. solutions in a 
weak sense, which are even more general than distributional solutions). An 
analogous result holds in the relativistic case, as was shown by Dudyhsky and 
Ekiel-Jezewska [SHI. Regarding classical solutions, Illner and Shinbrot |58| 
have shown global existence of solutions to the nonrelativistic Boltzmann 
equation for small initial data (close to vacuum). At present there is no 
analogous result for the relativistic Boltzmann equation and this must be re- 
garded as an interesting open problem. There is however a recent result [2] 
for the homogeneous relativistic Boltzmann equation where global existence 
is shown for small initial data and bounded scattering kernel. When the 
data are close to equilibrium (see below), global existence of classical solu- 
tions has been proved by Glassey and Strauss |Uj in the relativistic case 
and by Ukai |117j in the nonrelativistic case (see also |11U| and [69 ) . 

The collision operator Q(f,g) may be written in an obvious way as 

Q(f,g) = Q + (f,g)-Q~(f,g), 
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where Q + and Q~ are called the gain and loss term respectively. If the 
loss term is deleted the gain-term-only Boltzmann equation is obtained. It 
is interesting to note that the methods of proof for the small data results 
mentioned above concentrate on gain-term only equations and once that is 
solved it is easy to include the loss term. In [3J it is shown that the gain-term 
only classical and relativistic Boltzmann equations blow up for initial data 
not restricted to a small neihgbourhood of trivial data. Thus, if a global 
existence proof for unrestricted data will be given it will necessarily use the 
full collision operator. 

The gain term has a nice regularizing property in the momentum vari- 
able. In it is proved that given / G L 2 (R 3 ) and g £ L 1 ^ 3 ) with f,g>0, 
then 

\\Q + {f,g)\\m(Rf,) < C||/IIl2(r3 ) ||5'||li(r3 ) , (11) 

under some technical requirements on the scattering kernel. Here H s is the 
usual Sobolev space. This regularizing result was first proved by P.L. Li- 
ons in the classical situation. The proof relies on the theory of Fourier 
integral operators and on the method of stationary phase, and requires a 
careful analysis of the collision geometry, which is very different in the rel- 
ativistic case. See also |121j and |122j for a simplified proof in the classical 
and relativistic case respectively. 

The regularizing theorem has many applications. An important appli- 
cation is to prove that solutions tend to equilibrium for large times. More 
precisely, Lions used the regularizing theorem to prove that solutions to the 
(classical) Boltzmann equation, with periodic boundary conditions, converge 
in L 1 to a global Maxwellian, 

M = e -°lpl 2 +/3-P+7 a )7 Gfl, a>0,/3GM 3 , 

as time goes to infinity. This result had first been obtained by Arkeryd ^U] 
by using non-standard analysis. It should be pointed out that the conver- 
gence takes place through a sequence of times tending to infinity and it is 
not known whether the limit is unique or depends on the sequence. In the 
relativistic situation, the analogous question of convergence to a relativistic 
Maxwellian, or a Jiittner equilibrium solution, 

J = e -aVl+bl 2 +/3-p+ 7) Qj p and 7 ag abov6i with a > I/?), 

had been studied by Glassey and Strauss [37], In the periodic case 

they proved convergence in a variety of function spaces for initial data close 
to a Jiittner solution. Having obtained the regularizing theorem for the 
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relativistic gain term, it is a straightforward task to follow the method of 
Lions and prove convergence to a local Jiittner solution for arbitrary data 
(satisfying the natural bounds of finite energy and entropy) that is periodic 
in the space variables. In 11 it is next proved that the local Jiittner solution 
must be a global one, due to the periodicity of the solution. 

For more information on the relativistic Boltzmann equation on Minkowski 
space we refer to [101 > |49| , and |112j and in the nonrelativistic case we refer 
to the excellent review paper by Villani |118j and the books [IU], |24| . 

1.2 The Vlasov-Maxwell- and Vlasov-Poisson systems 

Let us consider a collisionless plasma, which is a collection of particles 
for which collisions are relatively rare and the interaction is through their 
charges. We assume below that the plasma consists only of one type of par- 
ticle, whereas the results below also hold for plasmas with several particle 
species. The particle rest mass is normalized to one. In the kinetic frame- 
work, the most general set of equations for modelling a collisionless plasma 
is the relativistic Vlasov-Maxwell system: 

dtf + v- V x f + (E(t, x)+vx B(t, x)) ■ V v f = (12) 



d t E + j = cV x B, V -E = p, (13) 
t B = -cV X E, V-B = 0. (14) 

The notation follows the one already introduced with the exception that the 
momenta are now denoted by v instead of p. This has become a standard 
notation in this field. E and B are the electric and magnetic fields, and v is 
the relativistic velocity 

v = — I =, (15) 

where c is the speed of light. The charge density p and current j are given 
by 

P= f dv , 3=1 vfdv. (16) 



Equation (|12|) is the relativistic Vlasov equation and ()13|) and (|14|) are the 
Maxwell equations. 

A special case in three dimensions is obtained by considering spherically 
symmetric initial data. For such data it can be shown that the solution also 
will be spherically symmetric, and that the magnetic field has to be constant. 
The Maxwell equation Vx£ = —dfB then implies that the electric field is 
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the gradient of a potential 0. Hence, in the spherically symmetric case the 
relativistic Vlasov-Maxwell system takes the form: 



d t f + v-V x f + (3E(t,x)-V v f = 
E = V(p, A(p = p. 



(17) 
(18) 



Here (3 = 1, and the constant magnetic field has been set to zero, since 
a constant field has no significance in this discussion. This system makes 
sense for any initial data, without symmetry constraints, and is called the 
relativistic Vlasov-Poisson equation. Another special case of interest is the 
classical limit, obtained by letting c — > oo in (|12j ) -p4j ) . yielding: 



where (3 = 1. See Schaeffer |106| for a rigorous derivation of this result. This 
is the (nonrelativistic) Vlasov-Poisson equation, and = 1 corresponds 
to repulsive forces (the plasma case). Taking (3 = — 1 means attractive 
forces and the Vlasov-Poisson equation is then a model for a Newtonian 
self-gravitating system. 

One of the fundamental problems in kinetic theory is to find out whether 
or not spontaneous shock formations will develop in a collisionless gas, i.e. 
whether solutions to any of the equations above will remain smooth for all 
time, given smooth initial data. 

If the initial data are small this problem has an affirmative solution in 
all cases considered above (see EH) E3> an d For initial data 

unrestricted in size the picture is more involved. In order to obtain smooth 
solutions globally in time, the main issue is to control the support of the 
momenta 

Q(t) := sup{|t>| : 3(s, x) G [0, t] x M 3 such that f(s, x, v) ^ 0}, (21) 

i.e. to bound Q(t) by a continuous function so that Q(t) will not blow 
up in finite time. That such a control is sufficient for obtaining global 
existence of smooth solutions follows from well-known results in the different 
cases (see EH1> |54| . |14| and |41|). For the full three-dimensional relativistic 
Vlasov-Maxwell system, this important problem of establishing whether or 
not solutions will remain smooth for all time is open. There has been an 
increased activity during the last years aiming at a solution of this problem. 
Two new methods of proofs of the same theorem as in |45j are given in |6U| 
and . A different sufficient criteria for global existence is given by Pallard 




(19) 
(20) 



E = V(p, A(p = p, 
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in [75] and he also shows a new bound for the electromagnetic field in terms 
of Q(t) in |77]. In two space and three momentum dimensions, Glassey and 
Schaeffer |42| have shown that Q(t) can be controlled, which thus yields 
global existence of smooth solutions in that case (see also [HI]). 

The relativistic and nonrelativistic Vlasov-Poisson equations are very 
similar in form. In particular, the equation for the field is identical in the 
two cases. However, the mathematical results concerning the two systems 
are very different. In the nonrelativistic case Batt Jl] gave an affirmative 
solution 1977 in the case of spherically symmetric data. Pfaffelmoser (HH] 
(see also Schaeffer |1U9| ') was the first one to give a proof for general smooth 
data. He obtained the bound 

Q(t) < C(l + t) (51+5)/n , 

where 5 > could be taken arbitrarily small. This bound was later improved 
by different authors. The sharpest bound valid for j3 = 1 and /3 = — 1 has 
been given by Horst [SU and reads 

Qit) < C(l + i)log(2 + i). 

In the case of repulsive forces {(3 = 1) Rein jSJ has found the sharpest 
estimate by using a new identity for the Vlasov-Poisson equation, discovered 
independently by Illner and Rein |S7] and by Perthame |79| . Rein's estimate 
reads 

Q(t) < C(l + t) 2 / 3 . 

Independently and about the same time as Pfaffelmoser gave his proof, Li- 
ons and Perthame jOE] used a different method for proving global existence. 
To some extent their method seems to be more generally applicable to at- 
tack problems similar to the Vlasov-Poisson equation but which are still 
quite different (see [2], |61j). On the other hand, their method does not give 
such strong growth estimates on Q(t) as described above. For the relativis- 
tic Vlasov-Poisson equation, Glassey and Schaeffer showed in the case 
(3 = 1 that if the data are spherically symmetric, Q(t) can be controlled, 
which is analogous to the result by Batt mentioned above. Also in the case 
of cylindrical symmetry they are able to control Q(t), see 01]. If /5 = — 1 
it was also shown in 41 that blow-up occurs in finite time for spherically 
symmetric data with negative total energy. This system, however, is un- 
physical in the sense that it is not a special case of the Einstein- Vlasov 
system. Quite surprisingly, for general smooth initial data none of the tech- 
niques discussed above for the nonrelativistic Vlasov-Poisson equation apply 
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in the relativistic case. This fact is annoying since it has been suggested 
that an understanding of this equation may be necessary for understand- 
ing the three-dimensional relativistic Vlasov-Maxwell equation. However, 
the relativistic Vlasov-Poisson equation lacks the Lorentz invariance; it is a 
hybrid of a classical Galilei invariant field equation and a relativistic trans- 
port equation (|17|) , Only for spherical symmetric data is the equation a 
fundamental physical equation. The classical Vlasov-Poisson equation is on 
the other hand Galilean invariant. In 3j a different equation for the field 
is introduced that is observer independent among Lorentz observers. By 
coupling this equation for the field to the relativistic Vlasov equation, the 
function Q(t) may be controlled as shown in 3j. This is an indication that 
the transformation properties are important in studying existence of smooth 
solutions (the situation is less subtle for weak solutions, where energy es- 
timates and averaging are the main tools, see jSH], JUS] and [57j). Hence, 
it is unclear whether or not the relativistic Vlasov-Poisson equation will 
play a central role in the understanding of the Lorentz invariant relativistic 
Vlasov-Maxwell equation. 

We refer to the book by Glassey BIT for more information on the rela- 
tivistic Vlasov-Maxwell system and the Vlasov-Poisson equation. 



1.3 The Nordstrom- Vlasov system 



Before turning to the main theme of this review, i.e. the Einstein- Vlasov 
system, we introduce a system which has many mathematical features in 
common with the Vlasov-Maxwell system and which recently has been math- 
ematically studied for the first time. It is based on an alternative theory of 
gravity which was given by Nordstrom |7U| in 1913. By coupling this model 
to a kinetic description of matter the Nordstrom- Vlasov system results. In 
Nordstrom gravity the scalar field <j) contains the gravitational effects as 
described below. The Nordstrom- Vlasov system reads 







A, 



dp 



(22) 



dtf + p- V x f ■ 
Here 



(dt<f> + p ■ V x ct>) P + (1 + M 2 r 1/2 V^j • V p f = 0. (23) 



P 



P 



VT+W 

denotes the relativistic velocity of a particle with momentum p. The mass 
of each particle, the gravitational constant and the speed of light are all 
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normalized to one. A solution (f , (j)) of this system is interpreted as follows: 
The spacetime is a Lorentzian manifold with a conformally flat metric which, 
in the coordinates (t, x), takes the form 

= e 2 *diag(-l, 1,1,1). 

The particle distribution / defined on the mass shell in this metric is given 
by 

f(t,x,p) = rtt,x,e*p). (24) 

The first mathematical study of this system was undertaken by Calogero 
in |2Uj where static solutions are analyzed and where also more details on 
the derivation of the system are given. Although the Nordstrom- Vlasov 
model of gravity does not describe physics correctly (in the classical limit 
the Vlasov-Poisson system of Newtonian gravity |2J is however obtained) 
it can nevertheless serve as a platform for developing new mathematical 
methods. The system has some common features with the Einstein- Vlasov 
system (see the next section) but seems in many respects to be less involved. 
The closest relationship from a mathematical point of view is the Vlasov- 
Maxwell system; this is evident in [22] where weak solutions are obtained 
in analogy with |35| . in [23] where a continuation criterion for existence of 
classical solutions is established in analogy with |45j (an alternative approach 
is given in [7S]) and in [SI] where global existence in the case of two space 
dimensions is proved in analogy with [IS] ■ The spherically symmetric case is 
studied in [H] and cannot be directly compared to the spherically symmetric 
Vlasov-Maxwell system (i.e. the Vlasov-Poisson system) since the hyperbolic 
nature of the equations is still present in the former system while it is lost 
in the latter case. In it is shown that classical solutions exist globally 
in time for compactly supported initial data under the additional condition 
that there is a lower bound on the modulus of the angular momentum of 
the initial particle system. It should be noted that this is not a smallness 
condition and that the result holds for arbitrary large initial data satisfying 
this hypothesis. 

1.4 The Einstein- Vlasov system 

In this section we will consider a self-gravitating collisionless gas and we will 
write down the Einstein- Vlasov system and describe its general mathemat- 
ical features. Our presentation follows to a large extent the one by Rendall 
in |100j . We also refer to Ehlers [HZ] and Stewart |lllj for more background 
on kinetic theory in general relativity. 
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Let M be a four-dimensional manifold and let g a t, be a metric with 
Lorentz signature (—,+,+,+) so that (M, g^) is a spacetime. We use the 
abstract index notation, which means that g a \, is a geometric object and not 
the components of a tensor. See |119j for a discussion on this notation. The 
metric is assumed to be time-orientable so that there is a distinction between 
future and past directed vectors. The worldline of a particle with non-zero 
rest mass m is a timelike curve and the unit future-directed tangent vector 
v a to this curve is the four-velocity of the particle. The four-momentum p a 
is given by mv a . We assume that all particles have equal rest mass m and we 
normalize so that m = 1. One can also consider massless particles but we 
will rarely discuss this case. The possible values of the four-momentum are 
all future-directed unit timelike vectors and they constitute a hypersurface P 
in the tangent bundle TM, which is called the mass shell. The distribution 
function / that we introduced in the previous sections is a non-negative 
function on P. Since we are considering a collisionless gas, the particles travel 
along geodesies in spacetime. The Vlasov equation is an equation for / that 
exactly expresses this fact. To get an explicit expression for this equation we 
introduce local coordinates on the mass shell. We choose local coordinates 
on M such that the hypersurfaces t = x° = constant are spacelike so that 
t is a time coordinate and x J , j = 1,2,3 are spatial coordinates (letters in 
the beginning of the alphabet always take values 0, 1, 2, 3 and letters in the 
middle take 1,2,3). A timelike vector is future directed if and only if its 
zero component is positive. Local coordinates on P can then be taken as 
x a together with the spatial components of the four-momentum p a in these 
coordinates. The Vlasov equation then reads 

d t f + ^d x J-±rl b p* P %f = 0. (25) 

Here a, b = 0, 1, 2, 3 and j = 1, 2, 3 and T J ab are the Christoffel symbols. It is 
understood that p° is expressed in terms of pi and the metric g a b using the 
relation g a bP a p b = — 1 (recall that m = 1). 

In a fixed spacetime the Vlasov equation is a linear hyperbolic equation 
for / and we can solve it by solving the characteristic system, 

dX'/ds = ^, (26) 

dP i /ds = -K b ^. (27) 
In terms of initial data /o the solution to the Vlasov equation can be written 
f(x a , P l ) = /o^O^VXP^xV)), (28) 
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where X l (s, x a ,p % ) and P l (s, x a ,p l ) solve (J3DJ and (JUJ) and where x a ,p l ) = 
x l and P l (t,x a ,p l ) =p l . 

In order to write down the Einstein- Vlasov system we need to define 
the energy-momentum tensor T ab in terms of / and g a b- In the coordinates 
(x a ,p a ) on P we define 



where as usual p a = g a bP b and \g\ denotes the absolute value of the determi- 
nant of g. Equation 1)25(1 together with Einstein's equations 



then form the Einstein- Vlasov system. Here G a b is the Einstein tensor, 
R a b the Ricci tensor, R is the scalar curvature and A is the cosmological 
constant. In most of this review we will assume that A = but section 2.3 
is devoted to the case of non-vanishing cosmological constant (where also 
the case of adding a scalar field is discussed). We now define the particle 
current density 



Using normal coordinates based at a given point and assuming that / is 
compactly supported it is not hard to see that T a b is divergence-free which 
is a necessary compatability condition since G a b is divergence-free by the 
Bianchi identities. A computation in normal coordinates also shows that 
A^ a is divergence-free, which expresses the fact that the number of particles 
is conserved. The definitions of T a b and N a immediately give us a num- 
ber of inequalities. If V a is a future directed timelike or null vector then 
we have N a V a < with equality if and only if / = at the given point. 
Hence N a is always future directed timelike if there are particles at that 
point. Moreover, if V a and W a are future directed timelike vectors then 
TabV a W b > 0, which is the dominant energy condition. If X a is a spacelike 
vector then T a bX a X b > 0. This is called the non-negative pressure condi- 
tion. These last two conditions together with the Einstein equations imply 
that R a bV a V b > for any timelike vector V a , which is the strong energy 
condition. That the energy conditions hold for Vlasov matter is one reason 
that the Vlasov equation defines a well-behaved matter model in general rel- 
ativity. Another reason is the well-posedness theorem by Choquet-Bruhat 
for the Einstein- Vlasov system that we will state below. Before stating that 
theorem we will first discuss the initial conditions imposed. 




G a b '■— Rab 



-^RQab = 8nT ab + Ag ab 
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The data in the Cauchy problem for the Einstein- Vlasov system consist of 
the induced Riemannian metric g%j on the initial hypersurface S, the second 
fundamental form kij of S and matter data /o- The relations between a 
given initial data set (gij,kij) on a three-dimensional manifold S and the 
metric g^j on the spacetime manifold is that there exists an embedding ip of 
S into the spacetime such that the induced metric and second fundamental 
form of ip(S) coincide with the result of transporting (gij,kij) with if). For 
the relation of the distribution functions / and /o we have to note that 
/ is defined on the mass shell. The initial condition imposed is that the 
restriction of / to the part of the mass shell over ip(S) should be equal to 
/o (V ;_1 ) d( , ip)~ 1 )o(fi where <j) sends each point of the mass shell over ip(S), to 
its orthogonal projection onto the tangent space to ip(S). An initial data set 
for the Einstein- Vlasov system must satisfy the constraint equations, which 
read 



Here p = T a i,n a n b and j a = —h ah T\ jC n c where n a is the future directed 
unit normal vector to the initial hypersurface and h ab = g ab + n a n b is the 
orthogonal projection onto the tangent space to the initial hypersurface. In 
terms of /o we can express p and j by (j a satisfies n a j a = so it can 
naturally be identified with a vector intrinsic to S): 



can now state the local existence theorem by Choquet-Bruhat |3L)| for the 
Einstein- Vlasov system. 

Theorem 1 Let S be a 3-dimensional manifold, gij a smooth Riemannian 
metric on S, a smooth symmetric tensor on S and fo a smooth non- 
negative function of compact support on the tangent bundle TS of S. Sup- 
pose that these objects satisfy the constraint equations \29\) and \3(J\) . Then 
there exists a smooth spacetime (M,g a i,), a smooth distribution function f 
on the mass shell of this spacetime, and a smooth embedding ip of S into 
M which induces the given initial data on S such that g a b and f satisfy the 
Einstein- Vlasov system and 4>(S) is a Cauchy surface. Moreover, given any 



R - hjk ij + (trk) 2 
V 4 fcf-Vi(trfc) 



167T/3, 

8vrj/. 



(29) 
(30) 




Here |^<7| is the determinant of the induced Riemannian metric on S. We 
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other spacetime (M',g' b ), distribution function f and embedding t/j' satis- 
fying these conditions, there exists a diffeomorphism x from an open neigh- 
bourhood of ip(S) in M to an open neighbourhood of ip'iS) * n M 1 which 
satisfies \ ip = $ an d carries g a b and f to g' ab and f respectively. 

In this context we also mention that local existence has been proved for the 
Yang-Mills- Vlasov system in |31j and that this problem for the Einstein- 
Maxwell-Boltzmann system is treated in This result is however not 

complete, the nonnegativity of / is left unanswered. The paper also raises 
the question on what hypothesis from a physical point of view are reasonable 
on the scattering kernel in curved spacetimes, an issue which does not seem 
fully understood. 

A main theme in the following sections is to discuss special cases for 
which the local existence theorem can be extended to a global one. There 
are interesting situations when this can be achieved, and such global exis- 
tence theorems are not known for Einstein's equations coupled to other forms 
of phenomenological matter models, i.e. fluid models (see, however, |27j ) . 
In this context it should be stressed that the results in the previous sections 
show that the mathematical understanding of kinetic equations on a flat 
background space is well-developed. On the other hand the mathematical 
understanding of fluid equations on a flat background space (also in the 
absence of a Newtonian gravitational field) is not satisfying. It would be 
desirable to have a better mathematical understanding of these equations 
in the absence of gravity before coupling them to Einstein's equations. This 
suggests that the Vlasov equation is natural as matter model in mathemat- 
ical general relativity. 

2 Global existence theorems for the Einstein-Vlasov 
system 

In general relativity two classes of initial data are distinguished. If an iso- 
lated matter distribution is studied, the data are called asymptotically flat. 
The initial hypersurface is topologically M 3 and (since far away from the 
body one expects spacetime to be approximately flat) appropriate fall off 
conditions are imposed. Roughly, a smooth data set (gij,kij, /o) on R 3 is 
said to be asymptotically flat if there exist global coordinates x l such that 
as \x\ tends to infinity the components gij in these coordinates tend to 
5ij, the components kij tend to zero, fo has compact support and certain 
weighted Sobolev norms of gij — Sij and kij are finite (see jlUU| ). The sym- 
metry classes that admit asymptotically flatness are few. The important 
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ones are spherically symmetric and axially symmetric spacetimes. One can 
also consider a case in which spacetime is asymptotically flat except in one 
direction, namely cylindrical spacetimes. A lot of work has been done on 
the spherically symmetric case and this will be reviewed below. In the case 
of cylindrical symmetry it has been shown that if singularities form then the 
first one must occur at the axis of symmetry |38| . 

Spacetimes that possess a compact Cauchy hypersurface are called cos- 
mological spacetimes and data are accordingly given on a compact 3-manifold. 
In this case the whole universe is modelled and not only an isolated body. 
In contrast to the asymptotically flat case, cosmological spacetimes admit a 
large number of symmetry classes. This gives one the possibility to study in- 
teresting special cases for which the difficulties of the full Einstein equations 
are strongly reduced. We will discuss below cases for which the spacetime 
is characterized by the dimension of its isometry group together with the 
dimension of the orbit of the isometry group. 



2.1 Spherically symmetric spacetimes 

The study of the global properties of solutions to the spherically symmetric 
Einstein- Vlasov system was initiated by Rein and Rendall in 1990. They 
chose to work in coordinates where the metric takes the form 

dg 2 = _ e 2»(t,r) dt 2 + e 2X(t,r) dr 2 + ^2 + ^2 

where t G R, r > 0, 8 G [0, 7r], tp G [0, 2n]. These are called Schwarzschild 
coordinates. Asymptotic flatness is expressed by the boundary conditions 

lim X(t,r) = lim fi(t,r) = 0, Mt > 0. 

r— >oo r— >oo 

A regular centre is also required and is guaranteed by the boundary condition 

A(i,0) = 0. 

With 

x = (r sin cp cos 6, r sin cj) sin 9, r cos (p) 
as spatial coordinate and 

v i =p > + e A - 1 — - — 

r r 

as momentum coordinates the Einstein- Vlasov system reads 
dtf + e^' A V ■ V x f - (A t — + e^VrV / n>F)- ■ V,/ = 0, (31) 

Vi + M 2 r r 



15 



e~ 2A (2rA r - 1) + 1 = 8irr 2 p, (32) 
e _2A (2r// r + 1) - 1 = 8vrrV (33) 

The matter quantities are denned by 



p{t,x) = / yJl + \v\ 2 f{t lX ,v) dv, (34) 



P(t,x) = [ (^Yf(t,x,v) y=t= . (35) 

Let us point out that this system is not the full Einstein- Vlasov system. The 
remaining field equations, however, can be derived from these equations. 
See [5H] for more details. Let the square of the angular momentum be 
denoted by L, i.e. 

T I 1 2 1 |2 I \2 

L := \x\ \v\ — (x ■ v) . 

A consequence of spherical symmetry is that angular momentum is con- 
served along the characteristics of ()31j) . Introducing the variable 

x ■ v 
w = , 

r 

the Vlasov equation for / = f(t,r,w,L) becomes 

d t f + e»- x ^d r f - (X t w + e^firE - e^ x ^)d w f = 0, (36) 
E r^E 

where 



E = E(r, w, L) = yjl + w 2 + L/r 2 . 
The matter terms take the form 



^ poo poo 

p(t,r) = — / Ef(t,r,w,L) dwdL, (37) 

r J-oo JO 

7T f°° f°° W 2 

P&r) = - I I —f(t,r,w,L)dwdL. (38) 



oo JO 



E 



Let us write down a couple of known facts about the system ()32|) . (|3"3"|) . (ffi)|) . H37|) . 
and (|38|) . A solution to the Vlasov equation can be written as: 

f(t, r, w, L) = / (i?(0, t, r, w, L), W(0, t, r, w, L), L), (39) 



16 



where R and W are solutions to the characteristic system, 

^ = (u-\)(s,R) W (AQ) 

ds E{R,W,LY K ' 



dW 
ds 



-X t (s, R)W - e ( ^ A)(s ' fl Vr(s, R)E(R, W, L) (41) 
L 



+e 



C*-A)(a,fl). 



R*E(R,W,Ly 



such that the trajectory (R(s, t, r, w, L), W(s, t, r, w, L), L) goes through the 
point (r, iv, L) when s = t. This representation shows that / is nonnegative 
for all t > and that / < ||/o||oo- There are two known conservation laws 
for the Einstein- Vlasov system: conservation of the number of particles, 

poo / roc roo \ 

4vr 2 / e A( *' r) I / / f(t,r,w,L)dLdw) dr, 

JO \J-oo Jo J 

and conservation of the ADM mass 

M := 4vr / r 2 p(t,r)dr. (42) 



J o 

Let us now review the global results concerning the Cauchy problem that 
have been proved for the spherically symmetric Einstein- Vlasov system. As 
initial data we take a spherically symmetric, nonnegative, and continuously 
differentiable function /o with compact support that satisfies 

pr poo poo 

Att 2 / / / Ef {r,w,L)dwdLdr < -. (43) 
Jo J-ooJo 2 

This condition guarantees that no trapped surfaces are present initially. 
In [HE] it is shown that for such an initial datum there exists a unique, con- 
tinuously differentiable solution / with /(0) = fo on some right maximal 
interval [0, T). If the solution blows up in finite time, i.e. if T < oo then p(t) 
becomes unbounded as t — > T. Moreover, a continuation criterion is shown 
that says that a local solution can be extended to a global one provided 
the u-support of / can be bounded on [0, T) (in [HHJ they chose to work 
in the momentum variable v rather than w,L.). This is analogous to the 
situation for the Vlasov-Maxwell system where the function Q(t) was intro- 
duced for the ^-support. A control of the ^-support immediately implies 
that p and p are bounded in view of (|34[) and (|35|) . In the Vlasov-Maxwell 
case the field equations have a regularizing effect in the sense that deriva- 
tives can be expressed through spatial integrals, and it follows [GStl] that 
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the derivatives of / also can be bounded if the f-support is bounded. For 
the Einstein- Vlasov system such a regularization is less clear, since e.g. fj, r 
depends on p in a pointwise manner. However, certain combinations of sec- 
ond and first order derivatives of the metric components can be expressed 
in terms of matter components only-without derivatives (a consequence of 
the geodesic deviation equation). This fact turns out to be sufficient for 
obtaining bounds also on the derivatives of / (see |88| for details). By con- 
sidering initial data sufficiently close to zero, Rein and Rendall show that 
the f-support is bounded on [0, T) and the continuation criterion then im- 
plies that T = co. It should be stressed that even for small data no global 
existence result like this one is known for any other phenomenological mat- 
ter model coupled to Einstein's equations. The resulting spacetime in [5H] 
is geodesically complete and the components of the energy momentum ten- 
sor as well as the metric quantities decay with certain algebraic rates in 
t. The mathematical method used by Rein and Rendall is inspired by the 
analogous small data result for the Vlasov-Poisson equation by Bardos and 
Degond |12j . This should not be too surprising since for small data the 
gravitational fields are expected to be small and a Newtonian spacetime 
should be a fair approximation. In this context we point out that in it 
is proved that the Vlasov-Poisson system is indeed the nonrelativistic limit 
of the spherically symmetric Einstein- Vlasov system, i.e. the limit when the 
speed of light c — > oo. (In (HZj this issue is studied in the asymptotically flat 
case without symmetry assumptions.) Finally, we mention that there is an 
analogous small data result using a maximal time coordinate jlOOj instead 
of a Schwarzschild time coordinate. In these coordinates trapped surfaces 
are allowed in contrast to the Schwarzschild coordinates. 

The case with general data is more subtle. Rendall has shown jlUlj 
that there exists data leading to singular spacetimes as a consequence of 
Penrose's singularity theorem. This raises the question of what we mean 
by global existence for such data. The Schwarzschild time coordinate is 
expected to avoid the singularity and by global existence we mean that so- 
lutions remain smooth as Schwarzschild time tends to infinity. Even though 
spacetime might be only partially covered in Schwarzschild coordinates, a 
global existence theorem for general data would nevertheless be very im- 
portant since it is likely that it would constitute a central step for proving 
weak cosmic censorship. Indeed, if this coordinate system can be shown 
to cover the domain of outer communications and if null infinity could be 
shown to be complete then weak cosmic censorship would follow. A partial 
investigation for general data in Schwarzschild coordinates was done in [?J2] 
and in maximal-isotropic coordinates in [100] . In Schwarzschild coordinates 
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it is shown that if singularities form in finite time the first one must be at 
the centre. More precisely, if f(t,r,w,L) = when r > e for some e > 0, 
and for all t, w and L, then the solution remains smooth for all time. This 
rules out singularities of the shell-crossing type, which can be an annoying 
problem for other matter models (i.e. dust). The main observation in |92| 
is a cancellation property in the term 

in the characteristic equation Q41|) . In 100 a similar result is shown but here 
also an assumption that one of the metric functions is bounded at the cen- 
tre is assumed. However, with this assumption the result follows in a more 
direct way and the analysis of the Vlasov equation is not necessary which in- 
dicates that such a result might be true more generally. Recently, Dafermos 
and Rendall [33] have shown a similar result for the Einstein- Vlasov system 
in double-null coordinates. The main motivation for studying the system in 
these coordinates has its origin from the method of proof of the cosmic cen- 
sorship conjecture for the Einstein-scalar field system by Christodoulou 26 . 
An essential part of his method is based on the understanding of the forma- 
tion of trapped surfaces ( [2H]). The presence of trapped surfaces (for the 
relevant initial data) is then crucial in proving the completeness of future 
null infinity in In [32] the relation between trapped surfaces and the 
completeness of null infinity was strengthened; a single trapped surface or 
marginally trapped surface in the maximal development implies complete- 
ness of null infinty. The theorem holds true for any spherically symmetric 
matter spacetime if the matter model is such that "first" singularities nec- 
essarily emanate from the center (where the notion of "first" is tied to the 
casual structure). The results in jU^l and in |100j are not sufficient for con- 
cluding that the hypothesis of the matter needed in the theorem in [32] is 
satisfied since they concern a portion of the maximal development covered by 
particular coordinates. Therefore, Dafermos and Rendall [33] choose double- 
null coordinates which cover the maximal development and they show that 
the mentioned hypothesis is satisfied for Vlasov matter. 

In 93 a numerical study was undertaken of the Einstein- Vlasov system 
in Schwarzschild coordinates. A numerical scheme originally used for the 
Vlasov-Poisson system was modified to the spherically symmetric Einstein- 
Vlasov system. It has been shown by Rein and Rodewis [M] that the numer- 
ical scheme has the desirable convergence properties. (In the Vlasov-Poisson 
case convergence was proved in |108j . See also |39|). The numerical experi- 
ments support the conjecture that solutions are singularity-free. This can be 
seen as evidence that weak cosmic censorship holds for collisionless matter. 
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It may even hold in a stronger sense than in the case of a massless scalar 
field (see |25ll26 _). There may be no naked singularities formed for any reg- 
ular initial data rather than just for generic data. This speculation is based 
on the fact that the naked singularities that occur in scalar field collapse 
appear to be associated with the existence of type II critical collapse, while 
Vlasov matter is of type I. This is indeed the primary goal of their numerical 
investigation-to analyze critical collapse and decide whether Vlasov matter 
is type I or type II. 

These different types of matter are defined as follows. Given small initial 
data no black holes are expected to form and matter will disperse (which 
has been proved for a scalar field |2Hj and for Vlasov matter [EE])- For large 
data, black holes will form and consequently there is a transition regime 
separating dispersion of matter and formation of black holes. If we introduce 
a parameter A on the initial data such that for small A dispersion occurs and 
for large A a black hole is formed, we get a critical value A c separating these 
regions. If we take A > A c and denote by rriB{A) the mass of the black hole, 
then if rriB(A) — > as A — > A c , we have type II matter, whereas for type I 
matter this limit is positive and there is a mass gap. For more information 
on critical collapse we refer to the review paper by Gundlach |5U| . 

For Vlasov matter there is an independent numerical simulation by 
Olabarrieta and Choptuik |75| (using a maximal time coordinate) and their 
conclusion agrees with the one in [22] ■ Critical collapse is related to self 
similar solutions; Martin-Garcia and Gundlach jSH] have presented a con- 
struction of such solutions for the massless Einstein- Vlasov system by using 
a method based partially on numerics. Since such solutions often are re- 
lated to naked singularities, it is important to note that their result is for 
the massless case (in which case there is no known analogous result to the 
small data theorem in jSH]) and their initial data are not in the class that 
we have described above. 

We end this section with a discussion of the spherically symmetric Einstein- 
Vlasov-Maxwell system, i.e. the case considered above with charged parti- 
cles. Whereas the constraint equations in the uncharged case do not involve 
any problems to solve once the distribution function is given (and satisfies 
l|43J) ) the charged case is more challenging. In [73] it is shown that solutions 
to the constraint equations do exist for the Einstein- Vlasov-Maxwell system. 
In |72| local existence is then shown together with a continuation criterion 
and in |JJ the regularity theorem in [22] is generalized to this case. 
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2.2 Cosmological solutions 

In cosmology the whole universe is modelled and the "particles" in the ki- 
netic description are galaxies or even clusters of galaxies. The main goal 
again is to determine the global properties of the solutions to the Einstein- 
Vlasov system. In order to do so, a global time coordinate t must be found 
(global existence) and the asymptotic behaviour of the solutions when t 
tends to its limiting values has to be analyzed. This might correspond to 
approaching a singularity (e.g. the big bang singularity) or to a phase of 
unending expansion. Since the general case is beyond the range of current 
mathematical techniques, all known results are for cases with symmetry (see 
however [S] where to some extent global properties are established in the case 
without symmetry). 

There are several existing results on global time coordinates for solu- 
tions of the Einstein- Vlasov system. In the spatially homogeneous case it is 
natural to choose a Gaussian time coordinate based on a homogeneous hy- 
persurface. The maximal range of a Gaussian time coordinate in a solution 
of the Einstein- Vlasov system evolving from homogenous data on a compact 
manifold was determined in j^Sj. The range is finite for models of Bianchi 
IX and Kantowski-Sachs types. It is finite in one time direction and infinite 
in the other for the other Bianchi types. The asymptotic behaviour of solu- 
tions in the spatially homogeneous case has been analyzed in |104j and 105 j. 
In |lU4j . the case of massless particles is considered, whereas the massive case 
is studied in |105| . Both the nature of the initial singularity and the phase 
of unlimited expansion are analyzed. The main concern is the behaviour of 
Bianchi models I, II, and III. The authors compare their solutions with the 
solutions to the corresponding perfect fluid models. A general conclusion is 
that the choice of matter model is very important since for all symmetry 
classes studied there are differences between the collisionless model and a 
perfect fluid model, both regarding the initial singularity and the expand- 
ing phase. The most striking example is for the Bianchi II models, where 
they find persistent oscillatory behaviour near the singularity, which is quite 
different from the known behaviour of Bianchi type II perfect fluid models. 
In |1U5| it also is shown that solutions for massive particles are asymptotic 
to solutions with massless particles near the initial singularity. For Bianchi I 
and II it also is proved that solutions with massive particles are asymptotic 
to dust solutions at late times. It is conjectured that the same holds true 
also for Bianchi III. This problem is then settled by Rendall in 103 . 

All other results presently available on the subject concern spacetimes 
that admit a group of isometries acting on two-dimensional spacelike or- 
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bits, at least after passing to a covering manifold. The group may be 
two-dimensional (local U(l) x U(l) or T 2 symmetry) or three-dimensional 
(spherical, plane, or hyperbolic symmetry). In all these cases, the quotient 
of spacetime by the symmetry group has the structure of a two-dimensional 
Lorentzian manifold Q. The orbits of the group action (or appropriate quo- 
tients in the case of a local symmetry) are called surfaces of symmetry. 
Thus, there is a one-to-one correspondence between surfaces of symmetry 
and points of Q. There is a major difference between the cases where the 
symmetry group is two- or three-dimensional. In the three-dimensional case 
no gravitational waves are admitted, in contrast to the two-dimensional 
case. In the former case, the field equations reduce to ODEs while in the 
latter their evolution part consists of nonlinear wave equations. Three types 
of time coordinates that have been studied in the inhomogeneous case are 
CMC, areal, and conformal coordinates. A CMC time coordinate, t, is 
one where each hypersurface of constant time has constant mean curvature 
(CMC) and on each hypersurface of this kind the value of t is the mean cur- 
vature of that slice. In the case of areal coordinates, the time coordinate is 
a function of the area of the surfaces of symmetry (e.g. proportional to the 
area or proportional to the square root of the area). In the case of conformal 
coordinates, the metric on the quotient manifold Q is conformally flat. 

Let us first consider spacetimes (M, g) admitting a three-dimensional 
group of isometries. The topology of M is assumed to be R x 5 1 x F, with F 
a compact two-dimensional manifold. The universal covering F of F induces 
a spacetime (M,g) by M = R x 5 1 x F and g = p*g, where p : M —> M 
is the canonical projection. A three-dimensional group G of isometries is 
assumed to act on (M,g). If F = S 2 and G = 50(3), then (M,g) is called 
spherically symmetric; if F = T 2 and G = E% (Euclidean group), then 
(M, g) is called plane symmetric; and if F has genus greater than one and 
the connected component of the symmetry group G of the hyperbolic plane 
H 2 acts isometrically on F = H 2 , then (M,g) is said to have hyperbolic 
symmetry. 

In the case of spherical symmetry the existence of one compact CMC 
hypersurface implies that the whole spacetime can be covered by a CMC 
time coordinate that takes all real values |95l I18j . The existence of one 
compact CMC hypersurface in this case was proved later by Henkel 52 using 
the concept of prescribed mean curvature (PMC) foliation. Accordingly this 
gives a complete picture in the spherically symmetric case regarding CMC 
foliations. In the case of areal coordinates, Rein |S2] has shown, under a size 
restriction on the initial data, that the past of an initial hypersurface can 
be covered. In the future direction it is shown that areal coordinates break 
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down in finite time. 

In the case of plane and hyperbolic symmetry, Rendall and Rein showed 
in |95| and |82| . respectively, that the existence results (for CMC time and 
areal time) in the past direction for spherical symmetry also hold for these 
symmetry classes. The global CMC foliation results to the past imply that 
the past singularity is a crushing singularity since the mean curvature blows 
up at the singularity. In addition, Rein also proved in his special case with 
small initial data that the Kretschmann curvature scalar blows up when the 
singularity is approached. Hence, the singularity is both a crushing and a 
curvature singularity in this case. In both of these works the question of 
global existence to the future was left open. This gap was closed in 0, and 
global existence to the future was established in both CMC and areal coor- 
dinates. The global existence result for CMC time is partly a consequence of 
the global existence theorem in areal coordinates, together with a theorem 
by Henkel [£2] that shows that there exists at least one hypersurface with 
(negative) constant mean curvature. Also, the past direction was analyzed 
in areal coordinates and global existence was shown without any smallness 
condition on the data. It is, however, not concluded if the past singularity 
in this more general case without the smallness condition on the data is a 
curvature singularity as well. The question whether the areal time coordi- 
nate, which is positive by definition, takes all values in the range (0, oo) or 
only in (Rq,oo) for some positive i?o is also left open. In the special case 
in [S2], it is indeed shown that Rq = 0, but there is an example for vacuum 
spacetimes in the more general case of U(l) x U(l) symmetry where i?o > 0. 
This question was resolved by Weaver |12L)j . She proves that if spacetime 
contains Vlasov matter (i.e. / ^ 0) then R = 0. Her result applies to a more 
general case which we now turn to. 

For spacetimes admitting a two-dimensional isometry group the first 
study was done by Rendall 96 in the case of local U(l)xU(l) symmetry (or 
local T 2 symmetry) . For a discussion of the topologies of these spacetimes we 
refer to the original paper. In the model case the spacetime is topologically 
of the form R x T 3 , and to simplify our discussion later on we write down 
the metric in areal coordinates for this type of spacetime: 

g = e 2 ^- u \-adt 2 + d6 2 ) + e- 2U t 2 [dy + Hd9 + Mdt] 2 

+e 2U [dx + Ady + (G + AH)dB + (L + AM)dt] 2 . (44) 

Here the metric coefficients r],U,a,A,H,L and M depend on t and 6 and 
6,x,y G S 1 . In [HI CMC coordinat es are considered rather than areal co- 
ordinates. The CMC- and the areal coordinate foliations are both geomet- 
rically based time foliations. The advantage with a CMC approach is that 
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the definition of a CMC hypersurface does not depend on any symmetry as- 
sumptions and it is possible that CMC foliations will exist for rather general 
spacetimes. The areal coordinate foliation, on the other hand, is adapted to 
the symmetry of spacetime but it has analytical advantages that we will see 
below. 

Under the hypothesis that there exists at least one CMC hypersurface, 
Rendall proves, without any smallness condition on the data, that the past of 
the given CMC hypersurface can be globally foliated by CMC hypersurfaces 
and that the mean curvature of these hypersurfaces blows up at the past 
singularity. Again, the future direction was left open. The result in [55] 
holds for Vlasov matter and for matter described by a wave map (which is 
not a phenomenological matter model). That the choice of matter model is 
important was shown by Rendall 99 who gives a non-global existence result 
for dust, which leads to examples of spacetimes that are not covered by 
a CMC foliation. 

There are several possible subcases to the U(l) x U(l) symmetry class. 
The plane case where the symmetry group is three-dimensional is one sub- 
case and the form of the metric in areal coordinates is obtained by letting 
A = G = H = L = M = and U = logi/2 in Another subcase, 

which still admits only two Killing fields (and which includes plane sym- 
metry as a special case), is Gowdy symmetry. It is obtained by letting 
G = H = L = M = in (|44|). In [3], the author considers Gowdy sym- 
metric spacetimes with Vlasov matter. It is proved that the entire maximal 
globally hyperbolic spacetime can be foliated by constant areal time slices 
for arbitrary (in size) initial data. The areal coordinates are used in a direct 
way for showing global existence to the future whereas the analysis for the 
past direction is carried out in conformal coordinates. These coordinates 
are not fixed to the geometry of spacetime and it is not clear that the entire 
past has been covered. A chain of geometrical arguments then shows that 
areal coordinates indeed cover the entire spacetime. This method was ap- 
plied to the problem on hyperbolic and plane symmetry in jjj . The method 
in [1] was in turn inspired by the work ^S] for vacuum spacetimes where the 
idea of using conformal coordinates in the past direction was introduced. As 
pointed out in [7j, the result by Henkel [53] guarantees the existence of one 
CMC hypersurface in the Gowdy case and, together with the global areal 
foliation in [3], it follows that Gowdy spacetimes with Vlasov matter can 
be globally covered by CMC hypersurfaces as well (also to the future). The 
general case of U(l) x U(l) symmetry was considered in [§j where it is shown 
that there exist global CMC and areal time foliations which complete the 
picture. In this result as well as in the preceeding subcases mentioned above 
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the question whether or not the areal time coordinate takes values in (0, oo) 
or in (R, oo), R > 0, was left open. This issue was solved by Weaver |12U| 
where she concludes that R = if the distribution function is not identically 
zero initially. 

A number of important questions remain open. To analyse the nature 
of the initial singularity, which at present is known only for small initial 
data in the case considered in {22^, would be very interesting. The question 
of the asymptotics in the future direction is also an important issue where 
very little is known. The only situation where a result has been obtained is 
in the case with hyperbolic symmetry. Under a certain size restriction on 
the initial data, Rein [HEj shows future geodesic completeness. However, in 
models with a positive cosmological constant more can be said. 

2.3 Cosmological models with a cosmological constant or a 
scalar field 

The present cosmological observations indicate that the expansion of the 
universe is accelerating and this has influenced the theoretical studies in 
the field during the last years. One way to produce models with accelerated 
expansion is to choose a positive cosmological constant in the Einstein equa- 
tions. Another way is to include a nonlinear scalar field among the matter 
fields. In this section we will review the recent results for the Einstein- Vlasov 
system where a cosmological constant or a linear- or nonlinear scalar field 
have been included into the model. 

As in the previous section we start with the models with highest degree 
of symmetry, i.e. the locally spatially homogeneous models. In the case of 
a positive cosmological constant Lee |62j has shown global existence as well 
as future causal geodesic completeness for initial data which have Bianchi 
symmetry. She also obtains the time decay of the components of the energy 
momentum tensor as t — > oo. The past direction for some spatially homoge- 
neous models is considered in |113j . Existence back to the initial singularity 
is proved and the case with a negative cosmological constant is discussed. 
In [S3] Lee considers the case with a nonlinear scalar field coupled to Vlasov 
matter. The form of the energy momentum then reads 



Here </> is the scalar field and V is a potential and the Bianchi identities lead 
to the following equation for the scalar field 



.Vlasov 



+ V^V^ - (-V^V 7 + V{(j)))g a p. 



(45) 



a/3 — 1 a/3 



v 7 </>v 7 </> = v'(4>). 



(46) 
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Under the assumption that V is nonnegative and C 2 global existence to the 
future is obtained and if the potential is restricted to the form 

V{<t>) = V e~ c *, 

where < c < 4-^/vr then future geodesic completeness is proved. 

In the previous section we discussed the situation when spacetime admits 
a three dimensional group of isometries and we distinguished three cases: 
plane, spherical and hyperbolic symmetry. In area time coordinates the 
metric takes the form 

ds 2 = _ e Kt,r) dt 2 + e X(t,r) dr 2 + ^2 + gj^^ 

where k = 0, —1, +1 correspond to the plane, spherical and hyperbolic case 
respectively and where sino# = l,shii# = sin# and sin_i# = sinh#. In |115j 
the Einstein- Vlasov system with a positive cosmological constant is investi- 
gated in the future (expanding) direction in the case of plane and hyperbolic 
symmetry. The authors prove global existence to the future in these coordi- 
nates and they also show future geodesic completeness. The positivity of the 
cosmological constant is crucial for the latter result. Recall that in the case 
of A = future geodesic completenss has only been established for hyper- 
bolic symmetry under a smallness condition of the initial data [S3] . Finally 
a form of the cosmic no hair conjecture is obtained for this class of space- 
times. Indeed, they show that the de Sitter solution acts as a model for the 
dynamics of the solutions by proving that the generalized Kasner exponents 
tend to 1/3 as t — > oo which in the plane case is the de Sitter solution. The 
remaining case of spherical symmetry is analysed in |114j . Recall that when 
A = 0, Rein [52] showed that solutions can only exist for finite time in the 
future direction in area time coordinates. By adding a positive cosmological 
constant global existence to the future is shown to hold true if initial data 
is given on t = to where i 2 , > 1/A. The asymptotic behaviour of the matter 
terms is also investigated and slightly stronger decay estimates are obtained 
in this case compared to the case of plane and hyperbolic symmetry. 

The results discussed so far in this section have concerned the future time 
direction and a positive cosmological constant. The past direction with a 
negative cosmological constant is analysed in |113| where it is shown that for 
plane and spherical symmetry the areal time coordinate takes all positive 
values which is in analogy with Weaver's |120j result for A = 0. If initial 
data is restricted by a smallness condition the theorem is proven also in the 
hyperbolic case and for such data the result of the theorem holds true in all 
of the three symmetry classes when the cosmological constant is positive. 
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The early-time asymptotics in the case of small initial data is also analysed 
and is shown to be Kasner-like. 

In |116j the Einstein- Vlasov system with a linear scalar field is analysed 
in the case of plane, spherical and hyperbolic symmetry. Here the potential 
V in (|45|) and (|46|) is zero. A local existence theorem and a continuation 
criterion, involving bounds on derivatives of the scalar field in addition to 
a bound on the support of one of the moment variables, is proven. For the 
Einstein-scalar field system, i.e. when / = 0, the continuation criterion is 
shown to be satisfied in the future direction and global existence follows in 
that case. 

3 Stationary solutions to the Einstein-Vlasov sys- 
tem 

Equilibrium states in galactic dynamics can be described as stationary so- 
lutions of the Einstein-Vlasov system, or of the Vlasov-Poisson system in 
the Newtonian case. Here, we will consider the former case for which only 
static, spherically symmetric solutions have been constructed, but we men- 
tion that in the latter case also, stationary axially symmetric solutions have 
been found by Rein [Mfij . 

In the static, spherically symmetric case, the problem can be formulated 
as follows. Let the space-time metric have the form 



where r > 0, £ [0, 7r], cp € [0, 2tt]. As before, asymptotic flatness is ex- 
pressed by the boundary conditions 

lim A(r) = lim /j,(r) = 0, Vt > 0, 

r—*oo r— »oo 

and a regular centre requires 

A(0) = 0. 

Following the notation in section 2.1, the time- independent Einstein-Vlasov 
system reads 



ds 2 



e Mr) dt 2 + e 2A(r) dr 2 + ^2 + ^2 



V 



= • V x / - VT+W e M_ V- • V,/ = 

2 V 



(47) 




(48) 
(49) 
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The matter quantities are denned as before: 



p(x) = / y 7 ! + \v\ 2 f(x,v) dv, (50) 




The quantities 

E := e fJl (r)y / l + \v\ 2 , L = \x\ 2 \v \ 2 - (x ■ v) 2 = \x x v\ 2 , 

are conserved along characteristics. E is the particle energy and L is the 
angular momentum squared. If we let 

f{x,v) = *{E,L), 

for some function the Vlasov equation is automatically satisfied. The form 
of $ is usually restricted to 

<S>(E,L) = </>(E)(L-L ) 1 , (52) 

where I > -1/2 and L > 0. If <p(E) = (E - E )%, k > -1, for some 
positive constant Eq this is called the polytropic ansatz. The case of isotropic 
pressure is obtained by letting / = so that <3? only depends on E. We refer 
to |HHj for information on the role of Lq. 

In passing, we mention that for the Vlasov-Poisson system it has been 
shown that every static spherically symmetric solution must have the 
form / = <&(E,L). This is referred to as Jeans' theorem. It was an open 
question for some time to decide whether or not this was also true for the 
Einstein- Vlasov system. This was settled in 1999 by Schaeffer jlU9j . who 
found solutions that do not have this particular form globally on phase 
space, and consequently, Jeans' theorem is not valid in the relativistic case. 
However, almost all results in this field rest on this ansatz. By inserting the 
ansatz for / in the matter quantities p and p, a nonlinear system for A and 
[i is obtained, in which 

e~ 2A (2rA r - 1) + 1 = %irr 2 Gq,{r, //), 
e~ 2X (2rfi r + 1) - 1 = 8vrr 2 i?$(r,/x), 

where 

27r r°° r 2 ^ 2 - 1 ) , , f 2 

G,(r,/i) = 4 He^e,L) dLde, 

r Ji Jo y'e^- 1 - L/r 

H*(r,fi) = -j / / $(ert r) e,L)y/e 2 - 1 - L/r 2 dLde. 

r Ji Jo 
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Existence of solutions to this system was first proved in the case of 
isotropic pressure in 90 and then extended to the general case in . 
The main problem is then to show that the resulting solutions have finite 
(ADM) mass and compact support. This is accomplished in .90 for a poly- 
tropic ansatz with isotropic pressure and in (SSj for a polytropic ansatz with 
possible anisotropic pressure. They use a perturbation argument based on 
the fact that the Vlasov-Poisson system is the limit of the Einstein- Vlasov 
system as the speed of light tends to infinity jH?5]. Two types of solutions are 
constructed, those with a regular centre ( [OH] and jHSj); and those ( 
with a Schwarzschild singularity in the centre. In 132] Rendall and Rein go 
beyond the polytropic ansatz and assume that $ satisfies 

§{E, L) = c(E - E ) k + L l + O((E - E) s + k )L l as E - E , 

where k > — 1, I > -1/2, k + / + 1/2 > 0, k < I + 3/2. They show that 
this assumption is sufficient for obtaining steady states with finite mass and 
compact support. The result is obtained in a more direct way and is not 
based on the perturbation argument mentioned above. Their method is 
inspired by a work on stellar models by Makino [HI], m which he considers 
steady states of the Euler-Einstein system. In |ST| there is also an interesting 
discussion about steady states that appear in the astrophysics literature. 
They show that their result applies to most of these steady states, which 
proves that they have the desirable property of finite mass and compact 
support. 

All solutions described so far have the property that the support of p 
contains a ball about the centre. In [HI] Rein shows that there exist steady 
states whose support is a finite, spherically symmetric shell, so that they 
have a vacuum region in the centre. 

At present, there are almost no known results concerning the stability 
properties of the steady states to the Einstein- Vlasov system. In the Vlasov- 
Poisson case, however, the nonlinear stability of stationary solutions has 
been investigated by Guo and Rein [HJ using the energy- Casimir method. 
In the Einstein- Vlasov case, Wolansky |123j has applied the energy-Casimir 
method and obtained some insights but the theory in this case is much 
less developed than in the Vlasov-Poisson case and the stability problem is 
essentially open. 
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